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Abstract. We provide here a simple, yet very general framework that allows 
us to explain several constraint propagation algorithms in a systematic way. In 
particular, using the notions commutativity and semi-commutativity, we show 
how the well-known AC-3, PC-2, DAC and DPC algorithms are instances of a 
single generic algorithm. The work reported here extends and simplifies that of 
Apti]. 



1 Introduction 

Constraint programming in a nutshell consists of formulating and solving so-called 
constraint satisfaction problems. One of the most important techniques developed in this 
area is constraint propagation that aims at reducing the search space while maintaining 
equivalence. 

We call the corresponding algorithms constraint propagation algorithms but sev- 
eral other names have also been used in the literature: consistency, local consistency, 
consistency enforcing. Waltz, filtering or narrowing algorithms. These algorithms usu- 
ally aim at reaching some form of "local consistency", a notion that in a loose sense 
approximates the notion of "global consistency". 

Over the last twenty few years several constraint propagation algorithms were pro- 
posed and many of them are built into the existing constraint programming systems. 
In Apt |jl|] we introduced a simple framework that allows us to explain many of these 
algorithms in a uniform way. In this framework the notion of chaotic iterations, so fair 
iterations of functions, on Cartesian products of specific partial orderings played a cru- 
cial role. In Monfroy and Rety [ pj[ | this framework was modified to study distributed 
chaotic iterations. This resulted in a general framework for distributed constraint prop- 
agation algorithms. 

We stated in Apt [|l]] that "the attempts of finding general principles behind the 
constraint propagation algorithms repeatedly reoccur in the literature on constraint sat- 
isfaction problems spanning the last twenty years" and devoted three pages to survey 
this work. Two references that are perhaps closest to our work are Benhamou and 
Telerman and Ushakov [|l6|. 

These developments led to an identification of a number of mathematical properties 
that are of relevance for the considered functions, namely monotonicity, inflationarity 
and idempotence (see, e.g., Saraswat, Rinard and Panangaden [Op and Benhamou and 



Older [^). Here we show that also the notions of commutativity and so-called semi- 
commutativity are important. 

As in Apt [jl[], to explain the constraint propagation algorithms, we proceed here 
in two steps. First, we introduce a generic iteration algorithm on partial orderings and 
prove its correctness in an abstract setting. Then we instantiate this algorithm with spe- 
cific partial orderings and functions. The partial orderings will be related to the consid- 
ered variable domains and the assumed constraints, while the functions will be the ones 
that characterize considered notions of local consistency in terms of fixpoints. 

This presentation allows us to clarify which properties of the considered functions 
are responsible for specific properties of the corresponding algorithms. The resulting 
analysis is simpler than that of Apt because we concentrate here on constraint prop- 
agation algorithms that always terminate. This allows us to dispense with the notion of 
fairness. On the other hand, we can now prove stronger results by taking into account 
the commutativity and semi-commutativity information. 

This article is organized as follows. First, in Section |2[ drawing on the approach of 
Monfroy and Rety 1 13 1, we introduce a generic algorithm for the case when the partial 
ordering is not further analyzed. Next, in Section |^, we refine it for the case when the 
partial ordering is a Cartesian product of component partial orderings and in Section 
^ explain how the introduced notions should be related to the constraint satisfaction 
problems. 

In the next four sections we instantiate the algorithm of Section ^ or some of its re- 
finements to obtain specific constraint propagation algorithms. In particular, in Section 
Ijwe derive algorithms for arc consistency and hyper-arc consistency. These algorithms 
can be improved by taking into account information on commutativity. This is done in 
Section^ and yields the well-known AC -3 algorithm. Next, in Section ^ we derive an 
algorithm for path consistency and in Section ^ we improve it, again by using informa- 
tion on commutativity. This yields the PC- 2 algorithm. 

In Section ^ we clarify under what assumptions the generic algorithm of Section 
H can be simplified to a simple for loop statement. Then we instantiate this simplified 
algorithm to derive in Section 10 the DAC algorithm for directional arc consistency and 
in Section [ill the DPC algorithm for directional path consistency. Finally, in Section |2| 
we briefly discuss possible future work. 

So we deal here only with the classical algorithms that establish (directional) arc 
consistency and (directional) path consistency and that are more than twenty, respec- 
tively ten, years old. However, several more "modern" constraint propagation algo- 
rithms can also be explained in this framework. In particular, in Apt [|l], page 203] we 
derived from a generic algorithm a simple algorithm that achieves the notion of rela- 
tional consistency of Dechter and van Beek [^. In turn, we can use the framework of 
Section ^t o derive the adaptive consistency algorithm of Dechter and Pearl [|[|. Now, 
Dechter^ showed that this algorithm can be formulated in a very general framework 
of bucket elimination that in turn can be used to explain such well-known algorithms 
as directional resolution, Fourier-Motzkin elimination, Gaussian elimination, and also 
various algorithms that deal with belief networks. 

Due to lack of space we do not define here formally the considered local consistency 
notions and refer the interested reader instead to the original papers or to Tsang [p^. 



2 Generic Iteration Algorithms 



Our presentation is completely general. Consequently, we delay the discussion of con- 
straint satisfaction problems till Section ^ In what follows we shall rely on the follow- 
ing concepts. 

Definition 1. Consider a partial ordering {D, \—) with the least element -L and a finite 
set of functions F := {fi, . . ., fk} on D. 

- By an iteration of F we mean an infinite sequence of values do , di , . . . defined 
inductively by 

do ■■= -L, 

where each ij is an element of [l..k]. 

— We say that an increasing sequence do \Z di \Z d2. . . of elements from D eventually 
stabilizes at d if for some j > we have di = dft?r i > j. □ 

In what follows we shall consider iterations of functions that satisfy some specific 
properties. 

Definition 2. Consider a partial ordering {D, \Z) and a function f on D. 

- f is caZZed inflationary ifx C f{x) for all x. 

— f is called monotomc if X 'O y implies f{x) ^ f{y) for all x^y. □ 

The following simple observation clarifies the role of monotonicity. The subsequent 
result will clarify the role of inflationarity. 

Lemma 1 (Stabilization). Consider a partial ordering (D, □ ) with the least element 
_L and a finite set of monotonic functions F on D. 

Suppose that an iteration of F eventually stabilizes at a common fixpoint d of the 
functions from F. Then d is the least common fixed point of the functions from F. 

Proof. Consider a common fixpoint e of the functions from F. We prove that d C e. Let 
do, di, . . .he the iteration in question. For some j > Owe have di = d for i > j. 

It suffices to prove by induction on i that di C e. The claim obviously holds for 
i = since do = -L- Suppose it holds for some i > 0. We have di+i — fj{di) for some 

j e [l..k]. 

By the monotonicity of fj and the induction hypothesis we get fj{di) □ fj{e), so 
di+i E e since e is a fixpoint of fj. □ 

We fix now a partial ordering {D, \Z) with the least element ± and a set of functions 
F :— {/i, . . ., fk} on D. We are interested in computing the least common fixpoint of 
the functions from F. To this end we study the following algorithm that is inspired by 
a similar algorithm of Monfroy and Rety Jl3[ ] . 



Generic Iteration Algorithm (gi) 

d:= _L; 
G := F; 
while G / do 

choose g E G; 

G:=G-{g}; 

G := GUupdate{G,g,dy, 

d := 9{d) 

od 

where for all G, g, d the set of functions update{G, g, d) from F is such that 

A. {/ e - G I f{d) = dA f{g{d)) g{d)} C update{G, g, d), 

B. g{d) = d implies that update{G, g, d) = 0. 

Intuitively, assumption A states that update{G, g, d) at least contains all the func- 
tions from F — G for which d is a fixpoint but g{d) is not. The idea is that such func- 
tions are repeatedly added to the set G. In turn, assumption B states that no functions 
are added to G in case the value of d, did not change. 

An obvious example of an update function that satisfies assumptions A and B is 

update{G, g,d):={f€F-G\ f{d) = dA f{g{d)) ^ g{d)}. 

However, this choice of the update function is computationally expensive because for 
each function / in — G we would have to compute the values f{g{dj) and f{d). In 
practice, we are interested in some approximations of the above update function. We 
shall deal with this matter in the next section. 

We now prove correctness of this algorithm in the following sense. 

Theorem 1 (GI) . 

(i) Every terminating execution of the GI algorithm computes in d a common fixpoint 
of the functions from F. 

(ii) Suppose that all functions in F are monotonic. Then every terminating execution 
of the GI algorithm computes in d the least common fixpoint of the functions from 
F. 

(Hi) Suppose that all functions in F are infiationary and that {D, C.) is finite. Then 
every execution of the GI algorithm terminates. 

Proof. 

(i) Consider the predicate / defined by: 

I-.^Vf &F-G f{d) = d. 

Note that / is established by the assignment G := F. Moreover, it is easy to check that 
/ is preserved by each while loop iteration. Thus / is an invariant of the while loop of 
the algorithm. Hence upon its termination 

(G = 0) A / 



holds, that is 



eF f{d) = d. 



(m) This is a direct consequence of (i) and the Stabihzation Lemina[TJ 

(iii) Consider the lexicographic ordering of the partial orderings {D, □) and (TV, <), 
defined on the elements of _D x JVhy 

(di, ni) <iex (^2, iff di □ ^2 or {di = c?2 and ni < n2). 

We use here the inverse ordering □ defined by: di □ ^2 iff '^2 E di and ^ di- 

Given a finite set G we denote by cardG the number of its elements. By assumption 
all functions in F are inflationary so, by virtue of assumption B, with each while loop 
iteration of the modified algorithm the pair 

{d, card G) 

strictly decreases in this ordering <iex- But by assumption (D, C ) is finite, so (_D, □) 
is well-founded and consequently so is [D x A/", <iex)- This implies termination. □ 

In particular, we obtain the following conclusion. 

Corollary 1 (GI) . Suppose that {D, \Z) is a finite partial ordering with the least 
element -L. Let F be a finite set of monotonia and inflationary functions on D. Then 
every execution of the GI algorithm terminates and computes in d the least common 
flxpoint of the functions from F. □ 

In practice, we are not only interested that the update function is easy to compute 
but also that it generates small sets of functions. Therefore we show how the function 
update can be made smaller when some additional information about the functions in 
F is available. This will yield specialized versions of the GI algorithm. First we need 
the following simple concepts. 

Definition 3. Consider two functions f,g on a set D. 

— We say that f and g commute if f{g{x)) = g{ f{x)) for all x. 

- We call f idempotent if f{f{x)) — f{x) for all X. □ 

The following result holds. 
Theorem 2 (Update). 

(i) Ifupdate{G, g, d) satisfles assumptions A andR, then so does the function 

update{G, g,d) — {g \ g is idempotent}. 

(ii) Suppose that for each g ^ F the set of functions Comm{g) from F is such that 

- g ^ Gomm{g), 

— each element of Gomm{g) commutes with g. 



Ifupdate{G, g, d) satisfies assumptions A andB, then so does the function 



update{G , g , d) — Comm{g). 
Proof. It suffices to establish in each case assumption A. 

{i) Suppose that g is idempotent. Then any function / such that f(g{d)) ^ g{d) differs 
from g. 

(ii) Consider a function / from F—G such that f{d) = d and f{g{d)) ^ g{d). Suppose 
that / S Comm{g). Then f{g(d.)) = g( f{d)) = g{d) which is a contradiction. So / ^ 
Comm{g). Consequently, / G update{G, g, d) — Comm{g) by virtue of assumption 
A f or update{G,g,d). □ 

We conclude that given an instance of the GI algorithm that employs a specific 
update function, we can obtain other instances of it by using update functions modified 
as above. Note that both modifications are independent of each other and therefore can 
be applied together. In particular, when each function is idempotent and the function 
Comm. satisfied the assumptions of (ii), then if update{G, g, d) satisfies assumptions 
A and B, then so does the function update{G, g, d) — {Comm{g) U {g}). 

3 Compound Domains 

In the applications we study the iterations are carried out on a partial ordering that is 
a Cartesian product of the partial orderings. So assume now that the partial ordering 
{D, E ) is the Cartesian product of some partial orderings {Di, C j), for i e [l..n], 
each with the least element ±i. So D = Di x ■ ■ ■ x Dn. 

Further, we assume that each function from F depends from and affects only certain 
components of D. To be more precise we introduce a simple notation and terminology. 

Defuiition4. Consider a sequence of partial orderings {Di, C i), . . ., (£)„, E „). 

- By a scheme (on n) we mean a growing sequence of different elements from [l..n]. 

- Given a scheme s := ii, . . on n we denote by (Dg, C g) the Cartesian product 
of the partial orderings {Di. , C i. ), for j G [1../]. 

- Given a function f on Dg we say that f is with scheme s and say that / depends 
on i if i is an element of s. 

- Given an n-tuple d := di, . . .,dn from D and a scheme s := ii, . . .,ii on n we 
denote by d[s] the tuple di^, . . ., d^. In particular, for j G [l..n] d[j] is the j-th 
element ofd. □ 

Consider now a function / with scheme s. We extend it to a function /+ from D to 
D as follows. Take d & D. We set 

f+{d) := e 

where e[s] = f{d[s]) and e[n — s] = d[n — s], and where n — s is the scheme obtained 
by removing from 1, . . ., n the elements of s. We call /+ the canonic extension of / to 
the domain D. 



So . . ., d„) = (ei, . . ., e„) implies di = for any i not in the scheme s of /. 

Informally, we can summarize it by saying that /+ does not change the components on 
which it does not depend. This is what we meant above by stating that each considered 
function affects only certain components of D. 

We now say that two functions, / with scheme s and g with scheme t commute if 
the functions /+ and commute. 

Instead of defining iterations for the case of the functions with schemes, we rather 
reduce the situation to the one studied in the previous section and consider, equivalently, 
the iterations of the canonic extensions of these functions to the common domain D. 
However, because of this specific form of the considered functions, we can use now a 
simple definition of the update function. More precisely, we have the following obser- 
vation. 

Note 1 ( Update ). Suppose that each function in F is of the form /+ . Then the following 
function update satisfies assumptions A and B: 

update{G, g^ ,d) 

{/+ £ F — G \ f depends on some i in s such that d[i] ^ g^{d)[i]}, 
where g is with scheme s. 

Proof. To deal with assumption A take a function /+ ^ F — G such that /+ {d) ~ d. 
Then /(e) — e for any e that coincides with d on all components that are in the scheme 
of/. 

Suppose now additionally that f^{g^{d)) ^ g^{d)- By the above g^{d) differs 
from d on some component i in the scheme of /. In other words, / depends on some i 
such that d[i] ^ g^{d)[i\. This i is then in the scheme of g. 

The proof for assumption B is immediate. □ 

This, together with the GI algorithm, yields the following algorithm in which we 
introduced a variable d' to hold the value of g+((i), and used F^ :— {/ | /+ G -F} and 
the functions with schemes instead of their canonic extensions to D. 

Generic Iteration Algorithm for Compound Domains (cd) 

d (J-i, . . ., ±„); 
d' := d; 
G i^o; 
while G 7^ do 

choose g E G; suppose g is with scheme s; 

G:=G-{5}; 

d'[s] -.^gidis])- 

G := G U {/ e Fq — G I / depends on some iin s such that d[i] ^ 
d[s] := d'[s] 

od 

The following corollary to the GI Theorem |l]and the Update Note |l] summarizes 
the correctness of this algorithm. 



Corollary 2 (CD) . Suppose that {D, \Z) is a finite partial ordering that is a Carte- 
sian product of n partial orderings, each with the least element _Lj with i G [l..n]. Let 
F be a finite set of functions on D, each of the form /+. 

Suppose that all functions in F are monotonic and infiationary. Then every execu- 
tion of the CD algorithm terminates and computes in d the least common fixpoint of the 
functions from F. □ 

In the subsequent presentation we shall deal with the following two modifications 
of the CD algorithm: 

- CD I algorithm. This is the version of the CD algorithm in which all the functions 
are idempotent and the function update defined in the Update Theorem ^(i) is 
used. 

- CDC algorithm. This is the version of the CD algorithm in which all the functions are 
idempotent and the combined effect of the functions update defined in the Update 
Theorem^ is used for some function Comm. 

For both algorithms the counterparts of the CD Corollary ^ hold. 



4 From Partial Orderings to Constraint Satisfaction Problems 

We have been so far completely general in our discussion. Recall that our aim is to 
derive various constraint propagation algorithms. To be able to apply the results of 
the previous section we need to relate various abstract notions that we used there to 
constraint satisfaction problems. 

This is perhaps the right place to recall the definition and to fix the notation. Con- 
sider a finite sequence of variables X := Xi,. . .,a;„, where n > 0, with respective 
domains T> :— Di, . . Dn associated with them. So each variable Xi ranges over the 
domain Di. By a constraint C on X we mean a subset of _Di x . . . x Dn. 

By a constraint satisfaction problem, in short CSP, we mean a finite sequence of 
variables X with respective domains V, together with a finite set C of constraints, each 
on a subsequence of X. We write it as (C ; xi S Di, . . .,Xn S Dn), where X := 
xi, . . .,Xn and 2? Di, . . .,Dn. 

Consider now an element d := di, . . ., c?„ of Di x ... x Dn and a subsequence 
Y := Xi-^, . . Xi^ of X. Then we denote by d\Y] the sequence di^, . . di^ . 

By a solution to (C ; xi e Di, . . .,a;„ e £)„) we mean an element d G Di x 
... X Dn such that for each constraint C G C on a sequence of variables Y we have 
d\Y] G C. We call a CSP consistent if it has a solution. Two CSP's Vi and V2 with the 
same sequence of variables are called equivalent if they have the same set of solutions. 
This definition extends in an obvious way to the case of two CSP's with the same sets 
of variables. 

Let us return now to the framework of the previous section. It involved: 

(i) Partial orderings with the least elements; 

These will correspond to partial orderings on the CSP's. In each of them the original 
CSP will be the least element and the partial ordering will be determined by the 
local consistency notion we wish to achieve. 



(ii) Monotonic and inflationary functions with schemes; 

These will correspond to the functions that transform the variable domains or the 
constraints. Each function will be associated with one or more constraints. 

(iii) Common fixpoints; 

These will correspond to the CSP's that satisfy the considered notion of local con- 
sistency. 

In what follows we shall discuss two specific partial orderings on the CSP's. In each 
of them the considered CSP's will be defined on the same sequences of variables. 

We begin by fixing for each set D a collection J-{D) of the subsets of D that in- 
cludes D itself. So is a function that given a set D yields a set of its subsets to which 
D belongs. 

When dealing with the hyper-arc consistency J'{D) will be simply the set V{D) of 
all subsets of D but for specific domains only specific subsets of D will be chosen. For 
example, to deal with the the constraint propagation for the linear constraints on integer 
interval domains we need to choose for J-{D) the set of all subintervals of the original 
interval D. 

When dealing with the path consistency, for a constraint C the collection T{C) 
will be also the set 7'(C) of all subsets of C. However, in general other choices may 
be needed. For example, to deal with the cutting planes method, we need to limit our 
attention to the sets of integer solutions to finite sets of linear inequalities with integer 
coefficients (see Apt [jl], pages 193-194]). 

Next, given two CSP's, (f) := {C ; xi € Di, . . .,Xn 6 Dn) and -0 ■— {C ; xi G 
D'-y,...,Xn G -D'l), we write (f) i' iff 

- D'^e T{Di) (and hence D[ C D,) for i e [l..n], 

- the constraints in C are the restrictions of the constraints in C to the domains 
D[,...,D'^. 

So ■(/; if i/; can be obtained from by a domain reduction rule and the domains 
of belong to the appropriate collections of sets J-{D). 

Next, given two CSP's, := (Ci, . . ., Ck ; V£) and V' (C^ , . . ., ; V£), we 
write Ec '0 iff 

- C'i e T{Ci) (and hence C- C C,) for i e [l..fc]. 

In what follows we call the domain reduction ordering and Cc the constraint 
reduction ordering. To deal with the arc consistency, hyper-arc consistency and direc- 
tional arc consistency notions we shall use the domain reduction ordering, and to deal 
with path consistency and directional path consistency notions we shall use the con- 
straint reduction ordering. 

We consider each ordering with some fixed initial CSP V as the least element. In 
other words, each domain reduction ordering is of the form 

{{V'\V\=d V'},\=d) 
and each constraint reduction ordering is of the form 



Note that (C ; xi e D'^,...,Xn e D'^) (C ; xi G £>'/,... ,x„ e iff 
D D^'fori e [l..n]. 

This means that for V = {C \ xi G Di, . . ., x„ e Z?„) we can identify the domain 
reduction ordering {{V' \ V 'P'}, Qd) with the Cartesian product of the partial or- 
derings {!F{Di), D), where i e [I..71]. Additionally, each CSP in this domain reduction 
ordering is uniquely determined by its domains and by the initial V. 

Similarly, 

{C[,.. ., C, ; V£) (Cr, . . ., C^' ; V£) iff D C'l for i e [l..k]- 

This allows us for V = (Ci, . . ., ; VE) to identify the constraint reduction ordering 
{{V I V Cc 'P'}, !^c) with the Cartesian product of the partial orderings {!F{Ci), 3 
), where i G [l..k]. Also, each CSP in this constraint reduction ordering is uniquely 
determined by its constraints and by the initial V. 

In what follows instead of the domain reduction ordering and the constraint reduc- 
tion ordering we shall use the corresponding Cartesian products of the partial orderings. 
So in these compound orderings the sequences of the domains (respectively, of the con- 
straints) are ordered componentwise by the reversed subset ordering D. Further, in each 
component ordering {T{D), D) the set D is the least element. 

Consider now a function / on some Cartesian product J-{Ei ) x . . . x T{Em)- Note 
that / is inflationary w.rt. the componentwise ordering 3 if for all {Xi, . . G 
T{Ei) X ... X T{E^) we have Y, C for all i G [l..m], where /(Xi, . . X,„) = 
(Fi, . . ., Ym)- 

Also, / is monotonic w.r.t. the componentwise ordering 3 if for all (Xi, . . ., Xm), 
{X[, . . .,X!^) G J'iEi) X ... X T{Ern) such that X, C X^ for afl i G [l..m], the 
following holds: if 

/(Xi,. . = (Fi, . . and fiX[, . . .^X'J = (F/, . . .,y^), 

then Y, C F/ for afl ie[l..m]- 

In other words, / is monotonic w.rt. 3 iff it is monotonic w.rt. C. This reversal of 
the set inclusion of course does not hold for the inflationarity notion. 

5 A Hyper-arc Consistency Algorithm 

We begin by considering the notion of hyper-arc consistency of Mohr and Masini JT^ 
(we use here the terminology of Marriott and Stuckey [|lO|). The more known notion 
of arc consistency of Mackworth is obtained by restricting one's attention to binary 
constraints. 

To employ the CD I algorithm of Section]^ we now make specific choices involving 
the items (i), (ii) and (iii) of the previous section. 

Re: (i) Partial orderings with the least elements. 

As already mentioned in the previous section, for the function T we choose the 
powerset function V, so for each domain D we put J-{D) :— V{D). 

Given a CSP V with the sequence Di, . . ., £)„ of the domains we take the domain 
reduction ordering with V as its least element. As already noted we can identify this 



ordering with the Cartesian product of the partial orderings {V{Di), D), where i e 
[l..n]. The elements of this compound ordering are thus sequences {Xi, . . .,X„) of 
respective subsets of the domains Di, . . ., D„ ordered componentwise by the reversed 
subset ordering D. 

Re: (ii) Monotonic and inflationary functions with schemes. 

Given a constraint C on the variables yi, . . .,yk with respective domains Ei, . . ., E^, 
we abbreviate for each j e [l..fc] the set {d[j] \ d G C} to 11^(0). Thus IZ^j(C) con- 
sists of all j-th coordinates of the elements of C. Consequently, nj{C) is a subset of 
the domain Ej of the variable yj. 

We now introduce for each i G [1..A;] the following function ttj on ■p(£'i) x • • • x 
V{Ek): 

TTi{Xi, . . .,Xk) := (Xi, . . .,Xi-i,Xl,Xi^i, . . .,Xk) 

where 

X'i := JTi(Cn(Xi X ••• xXfc)). 

That is, X- = {d[i] \ d G Xi x ■ ■ ■ x X/. and d G C}. Each function tt^ is associated 
with a specific constraint C. Note that X- CXi, so each function m boils down to a 
projection on the i-th component. 

Re: (iii) Common fixpoints. 

Their use is clarified by the following lemma that also lists the relevant properties 
of the functions tTj. 

Lemma 2 (Hyper-arc Consistency). 

(i) A CSP {C ; xi £ Di, . . ., .t„ £ _D„) is hyper-arc consistent iff {Di, . . ., Z?„) is a 
common fixpoint of all functions tt^ associated with the constraints from C. 

( ii) Each projection function Wi associated with a constraint C is 

- inflationary w.r.t. the componentwise ordering D, 

- monotonic w.r.t. the componentwise ordering 3, 

- idempotent □ 

By taking into account only the binary constraints we obtain an analogous charac- 
terization of arc consistency. The functions tti and 7r2 can then be defined more directly 
as follows: 

MX,Y) := {X',Y), 
where X' := {a e X \3b €Y {a,b) G C}, and 

7r2{X,Y) {X,Y'), 

where Y' := {b GY \ 3a G X (a, 6) G C}. 

Fix now a CSP V. By instantiating the CD I algorithm with 

Fq := {/ I / is a tt^ function associated with a constraint of V} 

and with each ±i equal to Di we get the HYPER-ARC algorithm that enjoys following 
properties. 



Theorem 3 (HYPER-ARC Algorithm). Consider a CSP V := {C ; xi £ Di, . . ., x„ G 
Dn) where each Di is finite. 

The HYPER-ARC algorithm always terminates. Let V' be the CSP determined by 
V and the sequence of the domains D[, . . ., D!^ computed in d. Then 

(i) V' is the \—d-least CSP that is hyper-arc consistent, 

( ii) V' is equivalent to V. □ 

Due to the definition of the \—d ordering the item (i) can be rephrased as follows. 
Consider all hyper-arc consistent CSP's that are of the form (C ; xi G D[, . . ., Xn G 
D'^) where D'^ C Di for i G [l..n] and the constraints in C are the restrictions of the 
constraints in C to the domains D[, . . D'^. Then among these CSP's V' has the largest 
domains. 



6 An Improvement: the AC -3 Algorithm 

In this section we show how we can exploit an information about the commutativity of 
the TT, functions. Recall that in Section |^ we modified the notion of commutativity for 
the case of functions with schemes. We now need the following lemma. 

Lemma 3 (Commutativity). Consider a CSP and two constraints of it, C on the vari- 
ables yi, . . .,yk and E on the variables zi, . . Zf. 

(i) For i, j G [l..k] the functions tt^ and Hj of the constraint C commute. 

( ii) If the variables yi and Zj are identical then the functions -Ki of C and ttj of E 
commute. □ 

Fix now a CSP. We derive a modification of the HYPER-ARC algorithm by in- 
stantiating this time the CDC algorithm. As before we use the set of functions Fo := 
{/ I / is a TTi function associated with a constraint of V} and each _Li equal to Di. Ad- 
ditionally we employ the following function Comm, where tt^ is associated with a con- 
straint C: 

Comm{TTi) := {ttj \ i ^ j and vTj is associated with the constraint C} 
U {-Kj I TVj is associated with a constraint E and 

the i-th variable of C and the j-th variable of E coincide}. 

By virtue of the Commutativity Lemma ^ each set Comm{g) satisfies the assump- 
tions of the Update Theorem 

By limiting oneself to the set of functions vri and tt2 associated with the binary 
constraints, we obtain an analogous modification of the corresponding arc consistency 
algorithm. 

Using now the counterpart of the CD Corollary ^ for the CDC algorithm we conclude 
that the above algorithm enjoys the same properties as the HYPER-ARC algorithm, that 
is the counterpart of the HYPER-ARC Algorithm Theorem ^ holds. 

Let us clarify now the difference between this algorithm and the HYPER-ARC al- 
gorithm when both of them are limited to the binary constraints. 



Assume that the considered CSP is of the form (C ; VE). We reformulate the above 
algorithm as follows. Given a binary relation R, we put 

{(6, a) | (a, 6) £ R}. 

For Fq we now choose the set of the tti functions of the constraints or relations from 
the set 

5*0 := {C I C is a binary constraint from C} 
U {C^ I C is a binary constraint from C}. 

Finally, for each tti function of some C G 5*0 on x, ?/ we define 

Comm{TTi) :={/!/ is the tti function of C'^} 

U {/ I / is the TTi function of some E ^ So on x, z where z ^ y}. 

Assume now that 

for each pair of variables x,y at most one constraint exists onx,y. (1) 

Consider now the corresponding instance of the CDC algorithm. By incorporating 
into it the effect of the functions tti on the corresponding domains, we obtain the fol- 
lowing algorithm known as the AC -3 algorithm of Mackworth 

We assume here that VE := xi <E Di, . . .,Xn G 

AC-3 Algorithm 

5*0 := {C I C is a binary constraint from C} 

U {C^ I C is a binary constraint from C}; 
S :— Sq; 
whUe 7^ do 

choose C E S; suppose C is on Xi, xj; 
A ■.= {aeD,\3be Dj (a, b) e C}; 
if Di changed then 

S :— S L) {C e 5o I C" is on the variables y, Xi where y ^ Xj} 

fi; 

S ■.= S-{C} 

od 

It is useful to mention that the corresponding reformulation of the HYPER-ARC 
algorithm differs in the second assignment to S which is then 

5 S* U {C G 5*0 I C" is on the variables y, z where y is Xi or z is Xi}. 

So we "capitalized" here on the commutativity of the corresponding projection 
functions tti as follows. First, no constraint or relation on Xi, z for some z is added 
to S. Here we exploited part (m) of the Commutativity Lemma |[ 

Second, no constraint or relation on Xj , Xi is added to S. Here we exploited part (i) 
of the Commutativity Lemma ^ because by assumption ([l|) is the only constraint 
or relation on Xj , Xi and its tti function coincides with the tt2 function of C. 



In case the assumption (|I|) about the considered CSP is dropped, the resulting algo- 
rithm is somewhat less readable. However, once we use the following modified defini- 
tion of Comm{TTi): 

Comm{TTi) :={/!/ is the tti function of some E E Sq on x, z where z ^ y} 

we get an instance of the CDC algorithm which differs from the AC-3 algorithm in 
that the qualification "where y ^ Xj" is removed from the definition of the second 
assignment to the set S. 



7 A Path Consistency Algorithm 

The notion of path consistency was introduced in Montanari It is defined for spe- 
cial type of CSP's. For simplicity we ignore here unary constraints that are usually 
present when studying path consistency. 

Definition 5. We call a CSP normalized if it has only binary constraints and for each 
pairx, y of its variables exactly one constraint on them exists. We denote this constraint 
by C^^y. □ 

Every CSP with only unary and binary constraints is trivially equivalent to a normal- 
ized CSP. Consider now a normalized CSP V . Suppose that V = (Ci, . . ., Cfc ; T>£). 

We proceed now as in the case of hyper-arc consistency. First, we choose for the 
function T the powerset function. For the partial ordering we choose the constraint 
reduction ordering of Section ^ or rather its counterpart which is the Cartesian product 
of the partial orderings (7'(Ci), 3), where i £ [l..fc]. 

Second, we introduce appropriate monotonic and inflationary functions with schemes. 
To this end, given two binary relations R and S we define their composition • by 

R-S:= {(a, h) I 3c ((a, c) G i?, (c, h) £ S*)}. 

Note that if i? is a constraint on the variables x, y and S a constraint on the variables 
y, z, then i? • S* is a constraint on the variables x, z. 

Given a subsequence x, y, z of the variables of V we now introduce three functions 

onV{C^.y) -y^ViC^^,) -y^ViCy,,): 

f^jP,Q,R) (P',Q,i?), 

where P' :=PnQ- R^, 

fl,{P,Q,R) := (P,Q',i?), 

where Q' Q (1 P ■ R, and 

flAP,Q,R) {P,Q,R'), 

where R' ■.= RnP'^ ■ Q. 

Finally, we introduce common fixpoints of the above defined functions. To this end 
we need the following counterpart of the Hyper-arc Consistency Lemma 0. 



Lemma 4 (Path Consistency). 

(i) A normalized CSP (Ci, . . .,C'k ; I^S) is path consistent iff (Ci, . . ., Cfe) is a com- 
mon fixpoint of all functions (/J y)^, (/^ and (fyz)^ associated with the sub- 
sequences X, y, z of its variables. 

(ii) The functions f^ y, /| ^ cind fy .^ are 

- inflationary w.r.t. the componentwise ordering 

- monotonic w.r.t. the componentwise ordering D, 

- idempotent. □ 
We now instantiate the GDI algorithm with the set of functions 

Fq :— {/ \ x,y,z is a subsequence of the variables of V and / e {fx,y^ fx.z^ fy,z}}^ 

n :— k and each Li equal to Ci. 

Call the resulting algorithm the PATH algorithm. It enjoys the following properties. 

Theorem4 {VATiiA\gorii\ai\). Consider a nonnalized CSP P := (Ci,...,Cfc; VS). 
Assume that each constraint d is finite. 

The PATH algorithm always terminates. Let V' :— (C(, . . ., ; }, where the 
sequence of the constraints C(, . . ., is computed in d. Then 

(i) V' is the \—c-least CSP that is path consistent, 

( ii) V is equivalent to V. □ 

As in the case of the HYPER-ARC Algorithm Theorem |]the item (z) can be rephrased 
as follows. Consider all path consistent CSP's that are of the form (C(, • ■ •, ; P'S) 
where C'i C d for i e [l..fc]. Then among them V' has the largest constraints. 



8 An Improvement: the PC-2 Algorithm 

As in the case of the hyper-arc consistency we can improve the PATH algorithm by 
taking into account the commutativity information. 

Fix a normalized CSP V. We abbreviate the statement "2;, y is a subsequence of the 
variables of "P" lo x < y. We now have the following lemma. 

Lemma 5 (Commutativity). Suppose that x ~< y and let z^u be some variables ofV 
such that {u, z} fl {x, y} ~ 0. Then the functions y and f"y commute. □ 

In other words, two functions with the same pair of variables as a subscript commute. 

We now instantiate the CDC algorithm with the same set of functions Fq as in Sec- 
tion 0. Additionally, we use the function Comm defined as follows, where x < y and 
where z ^ {x, ?/}: 

Comm{f'^ y) = {f}^y I u {x,y,z}}. 

Thus for each function g the set Comm{g) contains precisely to — 3 elements, 
where to is the number of variables of the considered CSP. This quantifies the maximal 



"gain" obtained by using the commutativity information: at each "update" stage of the 
corresponding instance of the CDC algorithm we add up to to — 3 less elements than in 
the case of the corresponding instance of the CD I algorithm considered in the previous 
section. 

By virtue of the Commutativity Lemma ^ each set Comm{g) satisfies the assump- 
tions of the Update Theorem H(m). We conclude that the above instance of the CDC 
algorithm enjoys the same properties as the original PATH algorithm, that is the coun- 
terpart of the PATH Algorithm Theorem^ holds. To make this modification of the PATH 
algorithm easier to understand we proceed as follows. 

Each function of the form f^y where x < y and u {x,y} can be identified with 
the sequence u, y of the variables. (Note that the "relative" position of u w.r.t. x and 
y is not fixed, so x, u, y does not have to be a subsequence of the variables of V.) This 
allows us to identify the set of functions Fq with the set 

Vo := {{x,u,y) \x ^y,u^ {x, y}}. 

Next, assuming that x y, we introduce the following set of triples of different 
variables of V: 

Vx,y ■■= {{x, y,u) \ X ^u}U {{y, x,u) \ y ^u} 
U {{u, x,y) \ u ^ y} U {{u, y,x) \ u ^ x}. 

Informally, Vx.y is the subset of Vb that consists of the triples that begin or end 
with either x, y or y, x. This corresponds to the set of functions in one of the following 
forms: f^y and fl^. 

The above instance of the CDC algorithm then becomes the following PC- 2 algo- 
rithm of Mackworth [^. Here initially E^.y = Cx,y 

PC-2 Algorithm 

Vo := {{x, u,y)\x <y,u(^ {x, y}}; 
V := Fo; 
while y 7^ do 

choose p eV\ suppose p = {x, u, y); 

apply y to its current domains; 

if Ex y changed tlien 
V:=VU Vx,y\ 

fi; 

V:=V- {p} 

od 

Here the phrase "apply /" y to its current domains" can be made more precise if the 
"relative" position of u w.r.t. x and y is known. Suppose for instance that u is "before" 
X and y. Then f^y is defined on V{Cu,x) x r{Cu,y) x V{Cx,y) by 

fx,yi-^i^,XT -^i^,tn -^x,y) i-^u,XT Eu,y, Ex,y n E^^^ ■ Eu,y), 



so the above phrase "apply f ^y to its current domains" can be replaced by the assign- 
ment 

.y • — . y ^ , X ' , y ■ 

Analogously for the other two possibilities. 

The difference between the PC-2 algorithm and the corresponding representation 
of the PATH algorithm lies in the way the modification of the set V is carried out. In 
the case of the PATH algorithm the second assignment to V is 

V:^VU V^.y U {{x, u, y) \ u ^ {x, y}}. 

9 Simple Iteration Algorithms 

Let us return now to the framework of Section ^. We analyze here when the while loop 
of the Generic Iteration Algorithm gi can be replaced by a for loop. First, we 
weaken the notion of commutativity as follows. 

Definition 6. Consider a partial ordering {D, \Z) and functions f and g on D. We say 
that f semi-commutes with g fw.rt. C ) if f{g{x)) C g{f{x)) for all x. □ 

The following lemma provides an answer to the question just posed. Here and else- 
where we omit brackets when writing repeated applications of functions to an argument. 

Lemma 6 (Simple Iteration). Consider a partial ordering (-D, C ) with the least ele- 
ment ±. Let F :— fi, . . fk be a finite sequence of monotonic, infiationary and idem- 
potent functions on D. Suppose that fi semi-commutes with fj for i > j, that is, 

h{fj{x))\=f,{f,{x))forallx. (2) 

Then /i/2. ■ .fki-L) is the least common fixpoint of the functions from F. □ 

Proof. We prove first that for i G [l..fc] we have 

f^flf2■..fki±)^flf2..■fk{l.). 

Indeed, by the assumption we have the following string of inclusions, where the last 
one is due to the idempotence of the considered functions: 

f^flf2■ ■ C /1/J2. . .fki±) E . . . C /1/2. . ./,/,. . □ /1/2. . 

Additionally, by the inflationarity of the considered functions, we also have for i G 
[l..k] 

/l/2-../fcW E/Jl/2.../fcW. 

So /1/2. . ./fc(^) is a common fixpoint of the functions from F. This means that 
the iteration of F that starts with ±, /fc(-L), fk-ifk{-L), . . ., /i/2- ■ -fki-L) eventually 
stabilizes at /i/2. ■ ./fc(^). By the Stabilization Lemma[l| we get the desired conclusion. 

□ 



The above lemma provides us with a simple way of computing the least common 
fixpoint of a set of finite functions that satisfy the assumptions of this lemma, in par- 
ticular condition (^. Namely, it suffices to order these functions in an appropriate way 
and then to apply each of them just once, starting with the argument _L. 

To this end we maintain the considered functions not in a set but in a list. Given a 
non-empty list L we denote its head by head(L) and its tail by tail(L). Next, given a 
sequence of elements ai, . . ., a„ with n > 0, we denote by [ai, . . ., a„] the list formed 
by them. If n = 0, then this list is empty and is denoted by [ ] and if n > 0, then 
head([ai, . . ., a„]) = ai and tail([ai, . . ., a„]) = [02, . . .a„]. 

The following algorithm is a counterpart of the GI algorithm. We assume in it that 
condition (Q) holds for the functions /i, . . ., /fc. 

Simple Iteration Algorithm (si) 
d := _L; 

L ■— [/fc, /fc-i, • ■ •, /i]; 
for i 1 to fc do 

g := head(i); 

L := tail(i); 

d g{d) 

od 

The following immediate consequence of the Simple Iteration Lemma ^ is a coun- 
terpart of the GI Corollary |l} 

Corollary 3 (SI ) . Suppose that (D, \Z) is a partial ordering with the least element 
_L. Let F := fi, . . fk be a finite sequence of monotonia, inflationary and idempotent 
fimctions on D such that holds. Then the S I algorithm terminates and computes in 
d the least common flxpoint of the functions from F. □ 

Note that in contrast to the GI Corollary |]we do not require here that the partial or- 
dering is finite. Because at each iteration of the for loop exactly one element is removed 
from the list L, at the end of this loop the list L is empty. Consequently, this algorithm 
is a reformulation of the one in which the line 

for I := 1 to fc do 

is replaced by 

while L ^ []do. 

So we can view the SI algorithm as a specialization of the GI algorithm of Section 
H in which the elements of the set of functions G (here represented by the list L) are 
selected in a specific way and in which the update function always yields the empty 
set. 

In Section |^ we refined the GI algorithm for the case of compound domains. An 
analogous refinement of the S I algorithm is straightforward and omitted. In the next 
two sections we show how we can use this refinement of the S I algorithm to derive two 
well-known constraint propagation algorithms. 



10 DAC: a Directional Arc Consistency Algorithm 



We consider here the notion of directional arc consistency of Dechter and Pearl [||]. To 
derive an algorithm that achieves this local consistency notion we first characterize it 
in terms of fixpoints. To this end, given a V and a linear ordering ^ on its variables, 
we rather reason in terms of the equivalent CSP obtained from V by reordering its 
variables along -< so that each constraint in is on a sequence of variables xi, . . .,Xk 
such that xi < X2 < ■ ■ ■ < Xk. 

The following characterization holds. 

Lemma 7 (Directional Arc Consistency). Consider a CSP V with a linear ordering 
-< on its variables. Let V-j (C ; xi G Di , . . .,Xn S Dn). Then V is directionally arc 
consistent w.rt. -< ijf{Di, . . ., Dn) is a common fixpoint of the functions tt^ associated 
with the binary constraints from "P^. □ 

We now instantiate in an appropriate way the S I algorithm for compound domains 
with all the tti functions associated with the binary constraints from V-^ . In this way 
we obtain an algorithm that achieves for V directional arc consistency w.r.t. First, 
we adjust the definition of semi-commutativity to functions with different schemes. To 
this end consider a sequence of partial orderings (Z^i, C j^), . . (Z?„, □ „) and their 
Cartesian product (Z), □ ). Take two functions, / with scheme s and g with scheme t. 
We say that / semi-commutes with g iw.r.t. □ ) if /+ semi-commutes with g+ w.r.t. 
C , that is if 

for all Q e D. 

The following lemma is crucial. 

Lemma 8 (Semi-commutativity). Consider a CSP and two binary constraints of it, 
Ci on M, z and C2 on x, y, where y ^ z. 

Then the tti function of Ci semi-commutes with the tti function of C2 w.r.t. the 
componentwise ordering D. □ 

Consider now a CSP V with a linear ordering -< on its variables and the corre- 
sponding CSP . To be able to apply the above lemma we order the tti functions of 
the binary constraints of in an appropriate way. Namely, given two vri functions, / 
associated with a constraint on u, z and g associated with a constraint on x, y, we put / 
before gify^z. 

More precisely, let x„ be the sequence of the variables of P^. So a;i -< X2 < 

. . . -< Xn- Let for m e [l..n] the list L„i consist of the tti functions of those binary 
constraints of that are on xj , x„i for some xj . We order each list i„, arbitrarily. 
Consider now the list L resulting from appending L„, L„-i, ■ ■ ., Li, in that order, so 
with the elements of i„ in front. Then by virtue of the Semi-commutativity Lemma ^ if 
the function / precedes the function g in the list L, then / semi-commutes with g w.r.t. 
the componentwise ordering D. 

We instantiate now the refinement of the S I algorithm for the compound domains 
by the above-defined list L and each _Li equal to the domain Di of the variable Xi. We 
assume that L has k elements. We obtain then the following algorithm. 



Directional Arc Consistency Algorithm (darc) 



d := (_Li, . . ., ±„); 
for j := 1 to A: do 

g := head(L); suppose g is with scheme s; 

L := tail(L); 

d[s] :=5(rfW) 

od 

This algorithm enjoys the following properties. 

Theorem 5 (DARC Algorithm). Consider a CSP V with a linear ordering -< on its 
variables. Let := (C ; xi G Di, . . ., G Dn)- 

The DARC algorithm always terminates. Let V be the CSP determined by and 
the sequence of the domains D'l, .... computed in d. Then 

(i) V' is the \—(i-least CSP in {Pi \ ^1} that is directionally arc consistent 
w.r.t. -<, 

( ii) V is equivalent to V. □ 

Note that in contrast to the HYPER- ARC Algorithm Theorem^ we do not need to 
assume here that each domain is finite. 

Assume now that for each pair of variables x, y of the original CSP V there exists 
precisely one constraint on x^y. The same holds then for V^. Suppose that := 
(C ; X\ € Di, . . .,Xn e Dn). Denote the unique constraint of on Xi,Xj by Cij. 
The above DARC algorithm can then be rewritten as the following algorithm known as 
the DAC algorithm of Dechter and Pearl 

for j := n to 2 by 1 do 
for i := 1 to J — 1 do 

A := {a G A I 3 6 e A (a- b) G C^^} 

od 

od 



11 DPC: a Directional Path Consistency Algorithm 

In this section we deal with the notion of directional path consistency defined in Dechter 
and Pearl [||]. As before we first characterize this local consistency notion in terms of 
fixpoints. To this end, as in the previous section, given a normalized CSP V we rather 
consider the equivalent CSP V^. The variables of are ordered according to -< and 
on each pair of its variables there exists a unique constraint. 

The following is a counterpart of the Directional Arc Consistency Lemma ^ 

Lemma 9 (Directional Path Consistency). Consider a normalized CSP V with a lin- 
ear ordering -< on its variables. Let :— (Ci, . . .,Ck "DS). Then V is directionally 
path consistent w.r.t. -< ijf (Ci, . . ., Cfe) is a common fixpoint of all functions (/J y)^ 
associated with the subsequences x, y, z of the variables ofV^. □ 



To obtain an algorithm that achieves directional path consistency we now instantiate 
in an appropriate way the S I algorithm. To this end we need the following lemma. 

Lemma 10 (Semi-commutativity). Consider a normalized CSP and two subsequences 
of its variables, Xi,yi, z and X2,y2, u. Suppose that u ^ z. 

Then the function f^^^y^ semi-commutes with the function fx^^y^ w.r.t. the compo- 
nentwise ordering D. □ 

Consider now a normalized CSP V with a linear ordering ~< on its variables and the 
corresponding CSP . To be able to apply the above lemma we order in an appropriate 
way the /* ^ functions, where the variables r, s, f are such that r -< s ^ t. Namely, we 
put /li.yi before f^^ y^ if u ^ z. 

More precisely, let xi, . . ., a;„ be the sequence of the variables of V-^, that is xi -< 
X2 < ■ ■ ■ -< Xn. Let for m £ [l..n] the list i,„ consist of the functions for 
some Xi and xj . We order each list L„i arbitrarily and consider the list L resulting from 
appending L„, L„_i, . . ., Li, in that order Then by virtue of the Semi-commutativity 
Lemma ^ if the function / precedes the function g in the list L, then / semi-commutes 
with g w.r.t. the componentwise ordering D. 

We instantiate now the refinement of the S I algorithm for the compound domains 
by the above-defined list L and each ±i equal to the constraint Ci. We assume that 
L has k elements. This yields the DIRECTIONAL Path Consistency Algorithm 
(DPATH) that, apart from of the different choice of the constituent partial orderings, 
is identical to the DIRECTIONAL ARC CONSISTENCY ALGORITHM DARC of the pre- 
vious section. Consequently, the DPATH algorithm enjoys analogous properties as the 
DARC algorithm. They are summarized in the following theorem. 

Theorem 6 (DPATH Algorithm). Consider a CSP V with a linear ordering -< on its 
variables. Let := (Ci, . . ., Cfe ; VS). 

The DPATH algorithm always terminates. Let V' :— (C(, . . ., ; "DE), where the 
sequence of the constraints C(, . . ., is computed in d. Then 

(i) V' is the \—c-least CSP in {Vi \ Vi} that is directionally path consistent 
w.r.t. 

( ii) V is equivalent to V. □ 

As in the case of the DARC Algorithm Theorem || we do not need to assume here 
that each domain is finite. 

Assume now that that xi , . . ., x„ is the sequence of the variables of V-^. Denote the 
unique constraint of V-i on Xi,Xj by Cij. 

The above DPATH algorithm can then be rewritten as the following algorithm known 
as the DPC algorithm of Dechter and Pearl [^: 

for TO := n to 3 by 1 do 
for j := 1 to TO — 1 do 
for i := 1 to J — 1 do 

od 

od 

od 



12 Conclusions 



In this article we introduced a general framework for constraint propagation. It allowed 
us to present and explain various constraint propagation algorithms in a uniform way. 
Using such a single framework we can easier verify, compare, modify, parallelize or 
combine these algorithms. The last point has already been made to large extent in Ben- 
hamou Additionally, we clarified the role played by the notions of commutativity 
and semi-commutativity. 

The line of research presented here could be extended in a number of ways. First, it 
would be interesting to find examples of existing constraint propagation algorithms that 
could be improved by using the notions of commutativity and semi-commutativity. 

Second, as already stated in Apt [jl]], it would be useful to explain in a similar way 
other constraint propagation algorithms such as the AC -4 algorithm of Mohr and Hen- 
derson [|ll|], the PC- 4 algorithm of Han and Lee [||], or the GAC-4 algorithm of Mohr 
and Masini [|l2|. The complication is that these algorithms operate on some extension 
of the original CSP. 

Finally, it would be useful to apply the approach of this paper to derive constraint 
propagation algorithms for the semiring-based constraint satisfaction framework of 
BistarelH, Montanari and Rossi that provides a unified model for several classes 
of "nonstandard" constraints satisfaction problems. 
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